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CN . Abstract 

I The theory of scale relativity provides a new insight into the origin 

" . of fundamental laws in physics. Its application to microphysics allows 

^ I ' us to recover quantum mechanics as mechanics on a non-differentiable 

, (fractal) spacetime. The Schrodinger and Klein-Gordon equations are 

■ demonstrated as geodesic equations in this framework. A development 

' of the intrinsic properties of this theory, using the mathematical tool of 

, Hamilton's bi-quaternions, leads us to a derivation of the Dirac equa- 

• tion within the scale-relativity paradigm. The complex form of the wave- 

I function in the Schrodinger and Klein-Gordon equations follows from the 

, non-differentiability of the geometry, since it involves a breaking of the 

^ ' invariance under the reflection symmetry on the (proper) time differential 

^ I ^ , element {ds <-> —ds). This mechanism is generalized for obtaining the 

f— I . bi-quaternionic nature of the Dirac spinor by adding a further symmetry 

' breaking due to non- different lability, namely the differential coordinate 

^ , reflection symmetry (da;^ <-> —dx^) and by requiring invariance under 

parity and time inversion. The Pauli equation is recovered as a non- 

^ ' relativistic-motion approximation of the Dirac equation. 
•i-H . 

H ■ 1 Introduction 

The theory of scale relativity consists of generalizing to scale transformations 
the principle of relativity, which has been applied by Einstein to motion laws. 
It is based on the giving up of the assumption of spacetime coordinate differen- 
tiability, which is usually retained as an implicit hypothesis in current physics. 
Even though this hypothesis can be considered as roughly valid in the classical 
domain, it is clearly broken by the quantum mechanical behavior. It has indeed 
been pointed out by Feynman (see, e.g., pj) that the typical paths of quantum 
mechanics are continuous but non-differentiable. 

In the present paper, after a reminder about the Schrodinger and Klein- 
Gordon equations, we apply the scale-relativistic approach to the understanding 
of the nature of bi-spinors and of the Dirac equation in a spacetime represen- 
tation. A step in this direction has been made by Gaveau et al T who have 



1 



generalized Nelson's stochastic mechanics 3 to the relativistic-motion case in 
(1+1) dimensions. However, an analytic continuation was needed to obtain the 
Dirac equation in such a framework, and, furthermore, stochastic mechanics is 
now known to be in contradiction with standard quantum mechanics as concerns 
multitime correlations in repeated measurements 0]. 

In the framework of a fractal spacetime viewed as a geometric analogue 
of quantum mechanics (S) |S| , Ord has developed a generalized version of the 
Feynman chessboard model which allows him to recover the Dirac equation in 
(3+1) dimensions, without analytical continuation |3|HliniEI|- We develop 
in the present paper an approach involving also a fractal spacetime, but in a 
different way: namely the fractality of spacetime is derived from the giving up 
of its differentiability, it is constrained by the principle of scale relativity and 
the Dirac equation is derived as an integral of the geodesic equation . This 
is not a stochastic approach (the use of statistically defined quantities is derived 
and we do not use Fokker-Planck equations), so that it does not come under 
the contradictions encountered by stochastic mechanics. 

One of the fundamental reasons for jumping to a scale-relativistic descrip- 
tion of nature is to increase the generality of the first principles retained to 
construct the theory. From a principle known as Occam's razor, the laws and 
structures of nature are the most general and simple laws and structures that 
are physically possible. It culminates nowadays in Einstein's "general" rela- 
tivity and his description of gravitation as a manifestation of the Riemannian 
geometry of spacetime, obtained by giving up the Euclidean flatness hypothesis. 
However, Einstein's principle of relativity is not yet fully general, since it ap- 
plies to coordinate transformations that are continuous and, at least two times, 
differentiable. The aim of the theory of scale relativity is to bring to light laws 
and structures that would be the manifestation of more general transformations, 
namely continuous ones, either differentiable or not. The standard "general rel- 
ativistic" theory will be recovered as the special differentiable case, while new 
effects would be expected from the non-differentiable part. 

In the present work, we focus our attention on the microphysical scale of 
motion in a non-differentiable spacetime, in the framework of "Galilean scale 
relativity" |12j . i.e. we consider fractal power law dilations with a constant 
fractal dimension. Therefore, we recover, from the first principles of this theory, 
(i) the four main evolution equations of standard quantum physics as geodesic 
equations on a fractal space/spacetime, namely the Schrodinger, (free) Klein- 
Gordon, Pauli and (free) Dirac equations, and (ii) the nature of the wavefunc- 
tions that appear in these equations, which are, respectively, complex, spinor 
and bi-spinor quantities. Indeed, as we shall see, the non-differentiability has 
two main consequences: (i) the fractality of the geometry that generates the 
purely quantum terms in the equations and (ii) discrete symmetry breakings at 
the infinitesimal level that lead to successive algebra doublings in the description 
of physical quantities. 

The derivation of the Schrodinger equation is given in Sec. |21 where we 
actualize the former works dealing with this issue, proposing a more accurate 
interpretation of the nature of the transition from the non-differentiable (fractal 
quantum scales) to the differentiable (classical scales) domain, and carefully 
justifying the different key choices made at each main step of the reasoning. An 
analogous updating is proposed in Sec. |21 for the free Klein-Gordon equation 
acting on a complex wavefunction. We obtain, in Sec. 01 the bi-quaternionic 
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form of the same Klein-Gordon equation, from which the Dirac equation then 
the Pauh equation naturally proceed, as shown in Sees. |5l and IHl The appendix 
is devoted to the presentation and justification of the mathematical tools we 
use. 



2 The Schrodinger equation revisited 

The Schrodinger equation is derived, in the framework of scale relativity, as a 
geodesic equation in a fractal three-space (see, e.g., (ElElEl)- In the present 
section, we update and give a more precise meaning to some of the most subtle 
issues we are dealing with. 

2.1 Transition from non-differentiability (fractal scales) to 
differentiability (classical scales) 

The aim of the present subsection is to identify the consequences of the giving 
up of the coordinate differentiability. We develop here the formalism relevant 
for the derivation of the non-relativistic-motion Schrodinger equation, valid in 
a fractal three-space with the time t as a curvilinear parameter. 

Strictly, the non-differentiability of the coordinates means that the velocity 

dt dt^a dt 

is undefined. This implies that, when dt tends to zero, either the ratio dX/dt 
tends to infinity, or it fluctuates without reaching any limit. This problem is 
overcome in the scale relativity framework thanks to the following fundamental 
theorem: a continuous and non-differentiable (or almost nowhere differentiable) 
curve is explicitly scale dependent, and its length tends to infinity when the 
scale interval tends to zero |12l 1181 ITHj . In other words, a continuous and non- 
differentiable space is fractal, in the general meaning given by Mandelbrot to 
this concept The velocity V{t), though it is undefined in the standard way, 
can now be defined in a new way as a "fractal function" ^(^j^^O; which 
is explicitly dependent on the scale interval 5t. As recalled in previous works 
|12l 1171 [T^ , these explicit scale variables can be defined only in a relative way 
(only their ratio has a physical meaning) so that the description of the scale 
space in which they lie is expected to come under a theory of scale relativity. 

In the case of a constant fractal dimension D, the relation between the 
position resolution interval and the time resolution interval is given by [121 116| 

5X w 5t^/". (2) 

The advantage of this method is that, for any given value of the resolution, 
5t, differentiability in t is recovered, which allows us to use the differential cal- 
culus, even when dealing with non-differentiability. However, the physical and 
the mathematical descriptions are not always coincident. Once 5t is given, one 
can write mathematical differential equations involving d/dt, make 9t — > 0, 
then solve for it and determine V(t^5t). Actually, this is a purely mathemati- 
cal description with no physical counterpart, since the very consideration of an 
interval dt < St, as occuring in an actual measurement, changes the function 
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V (such a behaviour, described by Heisenberg's uncertainty relations, is exper- 
imentally verified for any quantum system). However, making the particular 
choice dt = St induces a peculiar subspace of description where the physics 
and the mathematics coincide. We work, for the Schrodinger and Klein-Gordon 
equations, with such an identification of the time differential and of the new 
time resolution variable, but are led to give it up in Sees. 0] and [S] 

The scale dependence of the velocity suggests that we complete the standard 
equations of physics by new differential equations of scale. Writing the simplest 
possible equation for the variation of the velocity, V{t, dt), in terms of the new 
scale variable dt, as a first-order differential equation, dV/ din dt = f3{V), then 
Taylor expanding it, using the fact that < 1 (in relativistic- motion units 
c = 1), we obtain the solution ^1] 



V 
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l-l/D' 



(3) 



The velocity is now made of two independent terms of different order, since 
v/w « dt^^^/^ is an infinitesimal. We call v the "classical part" of the velocity 
(see below the definition of the classical part operator C£) and w its explicitly 
scale-dependent "fractal part" . The transition scale r and the scaled fluctuation 
C are chosen such that Cl{0 = and C£{C'^) = 1. 

We recognize here the combination of typical fractal behaviour with fractal 
dimension D and of a breaking of the scale symmetry at the scale transition 
r. In what follows, r will be identified with the Einstein-de Broglie scale of 
the system (r = h/E), since V v when dt ^ t (classical behaviour), and 
V K w when dt ^ t (fractal behaviour). Recalling that D — 2 plays the role 
of a critical dimension |12[ I13j . we recover Feynman's result in the asymptotic 
scaling domain, w oc {dt/T)~^/'^ which allows us to identify the fractal 
domain with the quantum one. In this paper, we consider only the case of 
fractal dimension 2. 

The above description strictly applies to an individual fractal trajectory. 
Now, one of the geometric consequences of the non-differentiability and of the 
subsequent fractal character of space/spacetime itself (not only of the trajecto- 
ries) is that there is an infinity of fractal geodesies relating any couple of points 
of this space/spacetime ^JEI- We therefore suggest ^21 that the description 
of a quantum mechanical particle could be reduced to the geometric properties 
of the set of fractal geodesies that corresponds to a given state of this "parti- 
cle" . In such an interpretation, the "particle" is not identified with a point mass 
which would follow the geodesies, but we expect its internal properties, mass, 
spin and charges, to be defined as geometric properties of the fractal geodesies 
themselves. As a consequence, any measurement is interpreted as a sorting out 
(or selection) of the geodesies by the measuring device E| ■ 

The transition scale appearing in Eq. © yields two distinct types of be- 
haviour of the system (particle) depending on the resolution at which it is con- 
sidered. Equation Q multiplied by dt gives the elementary displacement, dX, 
of the system as a sum of two terms 

dX = dx + d^, (4) 

d^ representing the "fractal part" and dx, the "classical part", defined as 

dx = V dt, (5) 
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which becomes, for D — 2, 

d^^aV2Vdt^^^, (7) 

with 21? = To = rw^ and CZ(a) = 0,C£{a'^) = 1. The scale 21? = tq wiU be 
subsequently identified with the Compton scale, h/mc, i.e. it gives the mass of 
the particle up to fundamental constants (see the argument in what follows). 
We note, from Eqs. to 0, that dx scales as dt, while d^ scales as dt^ . 
Therefore, the behaviour of the system is dominated by the term in the non- 
differentiable "fractal" domain (below the transition scale), and by the dx one 
in the differentiable "classical" domain (above the transition scale). 

Now, the Schrodinger, Pauli, Klein-Gordon and Dirac equations give results 
applying to measurements performed on quantum objects, but achieved with 
classical devices, on the classical/quantum interface. The microphysical scale 
at which a physical system is considered induces the sorting out of a bundle 
of geodesies, corresponding to the scale of the system, while the measurement 
process implies a smoothing out of the geodesic bundle coupled to a transition 
from the "fractal" to the "classical" domain. We therefore define an operator 
C£{ ), which we apply to the fractal variables or functions each time we are 
drawn to the "classical" domain where the dt behaviour dominates. The effect 
of C£ is to extract, from the fractal variables or functions to which it is applied, 
the "classical part" , i.e. the part scaling as dt. This justifies our above writting 
C£{() = and C£{C,'^) = 1, as it is straightforward, from Eq. Q, that the terms 
involving C scale as dt^^''^ and those involving scale as dt. This also led us to 
state, for the a dimensionless coefficient in Eq. Q, C£{a) — and C£{a?) = 1. 
Note the improvement of our new definition in terms of the classical part with 
respect to the previous interpretation in terms of an averaging process |12| . 



2.2 Differential-time symmetry breaking 

Another consequence of the non-differentiable nature of space (spacetime) is the 
breaking of local differential (proper) time reflection invariance. The derivative 
with respect to the time i of a differentiable function / can be written twofold 

di^^.^fit+dt)-m^^.^^m-fit-dt) ^ 

dt dt^o dt dt^a dt 

The two definitions are equivalent in the differentiable case. In the non- 
differentiable situation, both definitions fail, since the limits are no longer de- 
fined. In the framework of scale relativity, the physics is related to the behaviour 
of the function during the "zoom" operation on the time resolution 5t, here iden- 
tified with the differential clement dt, which is now considered as an independent 
variable. Two functions and f'_ are therefore defined as explicit functions of 
the two variables t and dt 



f^it,dt) = l^^±^^^^pl^, (9) 



/i(MO = /M)^^i^^. (10) 
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We pass from one to the other by the transformation dt ^ —dt (local 
differential-time reflection invariance), which was an implicit discrete symmetry 
of differentiable physics. The non-differentiable geometry implies that this sym- 
metry is now broken. When applied to the space coordinates, these definitions 
yield two velocities that are fractal functions of the resolution, V+[x{t),t,dt] 
and V- [x{t),t, dt]. In order to go back to the "classical" domain and derive the 
"classical" velocities appearing in Eq. 0, we smooth out each fractal geodesic 
in the bundle selected by the zooming process with balls of radius larger than 
T. This amounts to carrying out a transition from the non-differentiable to 
the differentiable domain and leads to define two "classical" velocity fields now 
resolution independent: V+[x{t),t,dt > t] = C£{V+[x{t),t, dt]) = v+[x{t),t] and 
V-[x{t),t,dt > r] = Ce{V-[x{t),t,dt]) = v^[x(t),t]. After the transition, there 
is no reason for these two velocities to be equal. While, in standard mechanics, 
the concept of velocity was one-valued, we introduce two velocities instead of 
one, even when going back to the "classical" domain. This two-valuedness of 
the velocity vector finds its origin in a breaking of the discrete time reflection 
invariance symmetry (dt <-> —dt), which is itself a mathematical consequence of 
non-differentiability. Therefore, if one reverses the sign of the time differential 
element, w+ becomes u_. A natural solution to this problem is to consider both 
{dt > 0) and {dt < 0) processes on the same footing, and to combine them in a 
unique twin process in terms of which the invariance by reflection is recovered. 
The information needed to describe the system is doubled with respect to the 
standard description. A simple and natural way to account for this doubling 
is to use complex numbers and the complex product. It is the origin of the 
complex nature of the wavefunction of quantum mechanics (see the appendix) . 



2.3 Covariant derivative operator 

Finally, we describe, in the scaling domain, the elementary displacements for 
both processes, dX±, as the sum of a classical part, dx± = v± dt, and a fluctu- 
ation about this classical part, d(^±, which is of zero classical part, C£{d^±) = 0. 

dX+{t) = v+ dt + d^+{t), 

dX^{t) =v^ dt + dC^{t). (11) 

Two classical derivatives, d/dtj^ and d/dt^, are defined, using the classical 
part extraction procedure. Applied to the position vector, x, they yield two 
classical velocities 

. , . . (12) 

As regards the fluctuations, the generalization to three dimensions of the 
fractal behaviour of Eq. © is written (for D = 2) 

Ci{di±, di±j) =±2V dt I, J = X, y, z, (13) 

The classical part of every crossed product d^±i d^±j, with i 7^ j, is null. This 
is due to the fact that, even if each term in the product scales as dt^^^, each of 
them behaves as an independent fractal fluctuation around its own classical part. 
Therefore, when we smooth out the geodesic bundle during the transition from 
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the fractal to the classical domain, we apply a process which is mathematically 
(not physically) equivalent to a stochastic "Wiener" process, and also more 
general, since we do not need any Gaussian distribution assumption. Thus, we 
can apply to the classical part of the d£_{t) product the property of the product 
of two independent stochastic variables: i.e. the classical part of the product is 
the product of the classical parts, and therefore here zero. 

To recover local differential time reversibility in terms of a new complex 
process |12| . we combine the two derivatives to obtain a complex derivative 
operator (see the appendix for a more detailed justification) 



dt 

Applying this operator to the position vector yields a complex velocity 





( d 




i / d 


--) 


^1 


ydt+ 




~ 2 \dt^ 


dt_) 



(14) 



The minus sign in front of the imaginary term is chosen here in order to 
finally obtain the Schrodinger equation in terms of V'- The reverse choice would 
give the Schrodinger equation for the complex conjugate of the wavefunction -0^, 
and would be therefore physically equivalent. The real part, F, of the complex 
velocity, V, represents the standard classical velocity, while its imaginary part, 
?7, is a new quantity arising from non-differentiability. At the usual classical 
limit, u+ = u_ = w, so that V = v and U = 0. 

Contrary to what happens in the differentiable case, the total derivative 
with respect to time of a fractal function f{x{t), t) of integer fractal dimension 
contains finite terms up to higher order [H?] 

df ^df ^ df dx, ^ 1 d^f dx.dx, ^ 1 d^f dx,dx,dXk ^ 

dt dt dxi dt 2 dxidxj dt 6 dxidxjdxt dt 

In our case, a finite contribution only proceeds from terms of Z?-ordcr, while 
lesser-order terms yield an infinite contribution and higher-order ones are neg- 
ligible. For a fractal dimension D — 2, the total derivative is written 

df df ^,dX 1 d^f dXidX, 

— ^ — +\/f. \ i ^. (17) 

dt dt dt 2 dxidxj dt 

Usually the term dXidXj/dt is infinitesimal, but here its classical part re- 
duces to Ci{d£^i d^j)/dt, which is now finite. Therefore, thanks to Eq. (|13(l . the 
last term of the classical part of Eq. (|17|l amoimts to a Laplacian, and we obtain 

£.(|+.,^V±Z.a)/. (18, 

Substituting Eqs. 1)18(1 into Eq. H14|l . we finally get the expression for the 
complex time derivative operator 12 

^ ^+V.V - iVA . (19) 
dt dt ^ ' 

The passage from standard classical (everywhere differentiable) mechanics to 
the new non-differentiable theory is then implemented by replacing the standard 
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time derivative d/dt by the new complex operator d /dt |12| . In other words, 
this means that d /dt plays the role of a "covariant derivative operator" , i.e. 
a tool that implements the form invariance of the equations. Note that in 
this replacement, one should be cautious and take into account the fact that 
this operator is a linear combination of first- and second-order derivatives, in 
particular when applying it to products and composed functions (see |20| and 
the appendix about the Leibniz rule for this operator). 



2.4 Covariant mechanics induced by scale laws 

Let us now recall how one generalizes the standard classical mechanics using this 
covariance. We assume that the classical part of the system can be characterized 
by a Lagrange function C{x^ V, i), from which an action S is defined. 

S= [ C{x,V,t)dt. (20) 

The Lagrange function and the action are now complex and are obtained 
from the classical Lagrange function L{x^ dx/dt, t) and the classical action S by 
replacing d/dt by d/dt. The stationary action principle applied on this complex 
action yields generalized Euler-Lagrange equations (see the appendix), 

= ^ (21) 

which are the equations one would obtain from applying the covariant deriva- 
tive operator {d/dt d/dt) to the standard Euler-Lagrange equations. This 
demonstrates the self-consistency of the approach and vindicates the use of com- 
plex numbers. Other fundamental results of standard classical mechanics are 
also generalized in the same way. In particular, assuming homogeneity of space 
in the mean leads one to define a generalized complex momentum given by 

(22) 

If we now consider the action as a functional of the upper limit of integration 
in Eq. (|2()|l . the variation of the action from one trajectory to another nearby 
one yields a generalization of another well-known relation of mechanics |21| 

V = VS. (23) 



2.5 Generalized Newton-Schrodinger equation 

Let us now consider the general case when the structuring external field is a 
scalar potential, The Lagrange function of a closed system, L — ^mv^ — $, 
is generalized, in the classical domain, as C{x,V,t) = ^mV^ — The Euler- 
Lagrange equations keep the form of Newton's equation of dynamics 

TO— V = -V$, (24) 
dt 

which is now written in terms of complex variables and complex operators. 
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In the case when there is no external field, the covariance is explicit, since 
Eq. H24I) takes the form of the equation of inertial motion 

d:V/dt = 0. (25) 

In both cases, the complex momentum V reads 

V = mV, (26) 

so that, from Eq. H23(l . the complex velocity V appears as a gradient, namely 
the gradient of the complex action, 

V VS/m. (27) 

We now introduce a complex wavefunction Tp which is nothing but another 
expression for the complex action S, 

^^g»5/So_ (28) 

The factor Sq has the dimension of an action and must be introduced at 
least for dimensional reasons. The "0 function is therefore related to the complex 
velocity appearing in Eq. (|27|l as follows 

V = — V(lnV'). (29) 
m 

The fundamental equation of dynamics H24() in terms of the new quantity tp 
takes the form 

j5o-^(Vln?A) = V$. (30) 
at 

Now d and V do not commute. However, there is a particular choice of the 
arbitrary constant Sq for which (i{'Vln'ip)/dt is nevertheless a gradient. 

Replacing d /dt by its expression, given by Eq. H19() . and replacing V by its 
expression in Eq. (|29|l . we obtain 



V$ = iSo 



VlmP-il — (VlnV'.V)(VlnV') +2?A(VlnV') 
m 



Using the remarkable identities 



(31) 



(Vln/)2 + Aln/ = ^ , (32) 

VA = AV, (33) 

V(V/)2 = 2(V/.V)(V/), (34) 

we obtain 

V = 2(Vln7/'.V)(Vln7/') + A(Vlni/;). (35) 

We recognize, on the right-hand side of this equation, the two terms of 
Eq. H31() . respectively, in factor of So/m and V. Therefore, the choice 

So = 2mV (36) 
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allows us to simplify this right-hand side, which becomes a gradient. The wave- 
function in Eq. H28() is therefore defined as 

and is a solution of the fundamental equation of dynamics, Eq. (|24|l . written as 

4v = -2X'v(i-^ln7/'-t-I?^l =-V$/m. (38) 
dt y at tp j 

This equation can now be integrated in a universal way which yields 

P^AV^ + iP^V-^^-O, (39) 

up to an arbitrary phase factor which may be set to zero by a suitable choice of 
the phase of ip. 

We are now able to enlighten the meaning of the choice Sq = 2m'D. We have 
seen that it is only under this particular choice that the fundamental equation 
of dynamics can be integrated. If we do not make this choice, the "0 function 
is a solution of a third-order, non-linear equation such that no precise physical 
meaning can be given to it. We therefore claim that this choice has a profound 
physical significance, since the meaning of '0 is directly related to the fact that it 
is a solution of the Schrodinger equation. In this equation. So is nothing but the 
fundamental action constant ?i, while T) defines the fractal/non- fractal transition 
(i.e. the transition from explicit scale dependence to scale independence), which 
is based on the constant A — 2I?/c. Therefore, the relation So — 2mT> becomes 
a relation between mass and a length-scale, which can be written 

Ac = — . (40) 
mc 

We recognize here the definition of the Compton length. Its profound mean- 
ing - i.e. up to the fundamental constants h and c, that of inertial mass itself 
- is thus given, in our framework, by the transition scale from explicit scale 
dependence to scale independence. It will be completely enlightened in the 
relativistic-motion case. This length scale is to be understood as a structure of 
scale space, not of standard space. 

We recover, in this case, the standard form of Schrodinger's equation 

—AiP + ih—ij^^tp. (41) 
2m ot 

New insight about the statistical meaning of the wavefunction (Born pos- 
tulate) can also be gained from the very construction of the theory. Indeed, 
while in standard quantum mechanics the existence of a complex wavefunction 
is set as a founding axiom, in the scale relativity framework it is related to the 
twin velocity field of the infinite family of geodesies with which the "particle" 
is identified, through Eq. (|29|l . V = — 2zI?V(ln-0). This opens the possibility of 
getting a derivation of Born's postulate in this context. 

This question has already been considered by Hermann 22 , who obtained 
numerical solutions of the equation of motion H24|) in terms of a large number 
of explicit trajectories (in the case of a free particle in a box). He constructed a 
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probability density from these trajectories and recovered in this way solutions of 
the Schrodinger equation without writing it and without using a wavefunction. 

Let us indeed derive the Born postulate in the one-dimensional stationary 
case. Consider the velocity field given by the real part V of the complex velocity 
V. It has been defined in such a way that it is identified with the classical velocity 
at the classical limit. Since this is a purely geometrical "object", the density 
of the fluid corresponding to this velocity field is a probability density P which 
satisfies the continuity equation, dP/dt + div(Py) = 0. 

Let us now set p = ipip\ i.e. V = -/pexp{iS/2m'D). By replacing ij) by 
this expression in the Schrodinger equation . and by using the identity 
V = 2I5VS', its imaginary part now reads 

|^+div(py)=0. (42) 

Since P and p are both solutions of a continuity equation, one can easily prove 
that their ratio K — P/p is a solution of the equation dK/dt = dK/dt + 
V.WK = 0. In the stationary one-dimensional case it is solved as K ^constant, 
and therefore the squared modulus of the wavefunction is identified with the 
probability density after a proper normalization. The general time-dependent 
three-dimensional case will be considered in a forthcoming work. 

3 Klein-Gordon equation 

In Sec. 121 the Schrodinger equation has been derived as a geodesic equation in 
a fractal three-space, for non-relativistic motion, in the framework of Galilean 
scale relativity. In the rest of this paper we shall be concerned with relativistic 
motion in the same framework of Galilean scale relativity and shall derive the 
corresponding free-particle quantum mechanical equations (Klein-Gordon's and 
Dirac's) as geodesic equations in a four-dimensional fractal spacetime. 

3.1 Relativistic-motion covariant derivative operator 

Most elements of the approach summarized in Sec. |21 remain correct in the 
relativistic-motion case, with the time, t, replaced by the proper time, s, as the 
curvilinear parameter along the geodesies. Now, not only space, but the full 
spacetime continuum is considered to be non-differentiable, thus fractal. We 
consider a small increment dX^ of a non-differentiable four-coordinate along 
one of the geodesies of the fractal spacetime. We can, as above, decompose 
dX^ in terms of a classical part Ci{dX^) = dx^ = v^ds and a fluctuation re- 
spective to this classical part, d^'', such that Ci{d^^) = 0, by definition. The 
non-differentiable nature of spacetime implies the breaking of the reflection in- 
variance at the infinitesimal level. If one reverses the sign of the proper time 
differential element, the classical part of the velocity v+ becomes We there- 
fore consider again both the {ds > 0) and {ds < 0) processes on the same footing. 
Then the information needed to describe the system is doubled with respect to 
the standard differentiable description and can be once more accounted for by 
the use of complex numbers. The new complex process, as a whole, recovers 
again the fundamental property of microscopic reversibility. 
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The elementary displacement along a geodesic of fractal dimension D = 2, 
respectively, for the (+) and (— ) processes, is then 

rf^± =rf±a;^+rfC± =w±ds + u^%/2Pdsi/^ (43) 

with d±x^ = v'^ds, d^±. = Uj_\/2'Dds^/^ and Uj_, a dimcnsionlcss fluctuation. 
The length-scale 21? is introduced for dimensional purpose. We define the clas- 
sical derivatives, d/ds+ and d/ds-, using the classical part extraction procedure 
exposed in Sec. 12.31 as 



(44, 

. i5s^0± \ OS / 



d 

ds± 

Once applied to a;^, they yield two classical four- velocities 
d ,, , ^ „ d 



x^'{s)^v1 -f-x^'{s) = vt. (45) 



They can be combined to construct a complex derivative operator 

d _l f d d \ i / d d 
ds 2 V ds-L ds_ / 2 V ds+ ds_ 



(46) 



The application of this operator to the position vector yields a complex 
four-velocity 

V- = ^x- = - = - . (47) 

ds 2 2 ^ ' 

As regards the fluctuations, the generalization to four dimensions of the 
fractal behaviour described in Eq. 143|l gives 

Cl{d^1d£,l) = T2T>r]i"'ds. (48) 

As noted in Sec. 12.31 each term on the left-hand side product behaves as 
an independent fractal fluctuation around its own classical part, which justifies 
a mathematical treatment analogous to that of a stochastic Wiener process. 
We make in the present paper the choice of a (-1-,—,—,—) signature for the 
Minkowskian metric of the classical spacetime, t]^'^ . Now, a diffusion (Wiener) 
process makes sense only in i?^ where the "metric" 77'^" should be positive defi- 
nite, if one wants to interpret the continuity equation satisfied by the probability 
density as a Kolmogorov equation Several proposals have been made to 

solve this problem |24[ I25|. which are equivalent in the end, and amount to 
transforming a Laplacian operator in into a Dalembertian. Namely, the two 
differentials of a function f{x, s) may be written 

£^^[l+v^d,TVd^8,)f. (49) 

As we only consider s-stationary functions, not explicitly depending on the 
proper time s, the complex covariant derivative operator reduces to 



ds 



= (V^ + iI?9^)a^. (50) 



The plus sign in front of the Dalembertian comes from the choice of the 
metric signature. 
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3.2 Geodesic equation 

To write the equation of motion, we use a generalized equivalence principle 
(identical to a strong covariance principle). We therefore obtain a geodesic 
equation in terms of the covariant derivative 

(51) 

as 

We now introduce the complex action according to 

dS — dyS dx" — —mc dx" (52) 

The complex four-momentum can thus be written as 

Vy = mc Vy ^ -dyS. (53) 

Now, the complex action, 5, characterizes completely the dynamical state 
of the particle, and we can introduce a complex wavefunction 

^^g»s/5o_ (54) 
It is linked to the complex four- velocity by Eq. (|53(l . which gives 

V, = — a.lnV. (55) 
mc 

3.3 Free-particle Klein-Gordon equation 

Replacing, in Eq. H51|l . the covariant derivative by its expression given by 
Eq. 1)50(1 and the complex four- velocity by that of Eq. 1(55(1 , we obtain 

- -^9'' In i:dady In - 9^9^ In = 0. (56) 

m^c"^ mc 

The particular choice, — Ti — 2mcV^ analogous to the one discussed in 
Sec. 12.51 allows us to simplify the left-hand side of Eq. 1(56(1 . using the following 
identity (which generalizes its three-dimensional counterpart of Eq. 1(35(1 ') 



1 fdf^\ ( 1 



[^^] = [d^\nij + -d^]d^d''\ni^. (57) 



Dividing by the constant factor T>^ , we obtain the equation of motion of the 
free particle under the form 

a^'f^Uo. (58) 



Therefore, the Klein-Gordon equation (without electromagnetic field) 

TT) f 

d'^d^i, + — = 0, (59) 
n 

becomes an integral of motion of the free particle, provided the integration 
constant is chosen equal to a squared mass term, rr?(? jh^ . The quantum be- 
haviour described by this equation and the probabilistic interpretation given to 
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ip is here reduced to the description of a free faU in a fractal spacetime, in anal- 
ogy with Einstein's general relativity where a particle subjected to the effect of 
gravitation is described as being in free fall in a curved spacetime. 

Another important property is the way the characteristic length of the frac- 
tal (quantum) to differentiable (classical) transition (the Compton length of the 
particle in its rest frame) and the mass term in the Klein-Gordon equation ap- 
pear. Here, as in the non-relativistic-motion case (Schrodinger), we identify the 
constant 21? with the Compton length h/mc, provided So is the fundamental 
action constant h, in order to obtain the motion equation under the form of a 
vanishing four-gradient. As for the mass term in the final Klein-Gordon equa- 
tion, it appears as a mere integration constant. However, it becomes connected 
to the transition length in the rest frame when matching these results with the 
non-relativistic description. 

4 Bi-quaternionic Klein-Gordon equation 

It has been known for a long time that the Dirac equation proceeds from the 
Klein-Gordon equation when written in a quaternionic form |26[ 12 7| . How- 
ever, no physical reason was known for this important mathematical property. 
We propose in the current section to introduce naturally, as a consequence of 
the non-differentiable geometry, a bi-quaternionic covariant derivative operator, 
leading to the definition of a bi-quaternionic velocity and wavefunction, which 
we use to derive the Klein-Gordon equation in a bi-quaternionic form. We use 
the quaternionic formalism, as introduced by Hamilton |28| . and further devel- 
oped by Conway (231221 (see also Synge PHI and Scheffers p?T]). 

4.1 Symmetry breakings and bi-quaternionic covariant deriva- 
tive operator 

Most of the approach described in SecElremains applicable. However, the main 
new features obtained in the now studied case proceed from a deeper description 
of the scale formalism. Namely we now consider a more general case involving 
the subsequent breaking of the symmetries: 
ds —ds 
dxf^ ^ —dx^ 

Indeed, we have up to now considered only the effect of non-differentiability 
on the total derivative d/ds (that amounts to d/dt in the non-relativistic case 
leading to the Schrodinger equation). Now the velocity fields of the geodesic 
bundles are functions of the coordinates, so that we are led to analyse also the 
physical meaning of the partial derivatives d/dx (we use only one coordinate 
variable in order to simplify the writing) in the decomposition d/ds — d/ds + 
{dx/ds)d/dx. Strictly speaking, df /dx does not exist in the non-differentiable 
case. We are therefore once again led to introduce fractal functions f{x,5x), 
explicitly dependent on the coordinate resolution interval, whose derivative is 
undefined only at the unobservable limit Sx 0. As a consequence of the very 
construction of the derivative, which needs two points to be defined (instead 
of one for the position and time coordinates), there are two definitions of the 
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partial derivative of a fractal function instead of one, namely 



df _ f{x + dx, dx) ~ f{x, dx) ^^^^ 
dx+ dx ' 

df f{x,dx) — f{x — dx,dx) 



dx - dx 



(61) 



They are transformed one into the other under the reflection dx <-> —dx. 

The x^^ <-> —x'^ symmetry corresponds to the parity P and time-reversal T 
symmetries the breaking of which is already taken into account in the standard 
definition of Dirac spinors in terms of pairs of Pauli spinors 32. . 

We have already stressed, in Sec. 12.11 that differentiability can be recovered 
(at large scales), even when dealing with non-differentiability. However, this 
implies that, for any set of differential equations describing a given process, the 
physical and mathematical descriptions are only coincident on a limited scale 
range. We can say that any consistent mathematical tool lives in a description 
space which is tangent to the physical space and that the validity of this tool 
is therefore limited to a finite scale region around the contact point. For the 
Schrodinger and complex Klein-Gordon equations, the tangent mathematical 
space is such that dt = St (Schrodinger) or ds = 6s (Klein-Gordon). When 
jumping to smaller (higher energy) scales, we are led to give up this peculiar 
choice, and retain, at the Dirac scale, a new mathematical description where the 
differentials and the resolution variables do no longer coincide. However, this 
description is only valid at the scales where the Dirac equation applies without 
corrections and we shall be led to improve it when going to yet smaller scales. 

In the scaling domain, the four spacetime coordinates X'^(s, e^, e^) are fractal 
functions of the proper time, s, and of the resolutions, and e^. We consider 
the case when, for an elementary displacement dX^^ corresponding to a shift ds 
in the curvilinear parameter, the resolutions verify < dX^^ and < ds. 

We can apply, to these two different elementary displacements, the reasoning 
developed in Sees. 12.21 and 12.31 which yields, considering first the displacement 
dX^, the canonical decomposition 

dX^' ^dx^" +d^'', (62) 

with 

CiidX") = dxf" = ds, (63) 
d^" = a^^V2P(ds2)A, Cl{a'-l) = 0, C£{{a1)^) = 1, (64) 

and 

- dXf" = fe" + 5^", (65) 

with 

C£{-dX'') = Sx^" = vt ds, (66) 

S^" = a''_V2D{ds^)^ , Ci{at) = 0, Ce{{a'lf) = 1. (67) 
In the differentiable case, dX^ = —{—dX'^), and therefore u+ — —vt . This 

is no longer the case in the non-differentiable case, where the local symmetry 
dx'^ <-» —dx^ is broken. 
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Furthermore, we must also consider the breaking of the symmetry ds ^ —ds 
proceeding from the twofold definition of the derivative with respect to the 
curvilinear parameter s. Applied to X^, considering an elementary displacement 
ds, the classical part extraction process gives two classical derivatives d/ds+ 
and d/ds^, which yield in turn two classical velocities, which we denote by 
w i- + . Considering now the same extraction process, applied to an elementary 

displacement —ds, the classical derivatives d/ ds+ and d/ds^ allow us again to 
define two classical velocities, denoted by w j_ _ . We summarize this result as 



<+=^, (68) 
s M ds± s p. ds± 

We can, at this stage, define several total derivatives with respect to s of 
a fractal function /. We write them - using a compact straightforward no- 
tation with summation over repeated indices - after substituting, in the four- 
dimensional analogue of Eq. H17() . the expressions for the derivatives of which 
are obtained when using Eqs. (|62|) to l|67|) with the expressions of Eq. (|68|l for 
the classical velocities 



-_±±±±± = +x;^^+u;'i± — — + Q^Q^P 

s X » ^ t ()s s M a M oX^^ dX^^dX'^ 



1 1 ^ 1 1 = + ± ± + < ± + "i^^'^^TTTi^ ' (69) 

with 

„2 



= V2X>(ds^)ro- 2 . (70) 

Now, when we apply the classical part operator to Eq. H69|) . using Eqs. (|64|l . 
(|67|l and (|7()|l . the w'^ disappear at the first order, but, at the second order, for 
the fractal dimension D = 2, we obtain the analogue of Eq. I|48(l 

Cl{w'i ±w'i±) = T'2'DT]'"'ds, (71) 

S fl Bp 

The =F sign on the right-hand side is the inverse of the s-sign on the left-hand 
side. We can thus write the classical part of the total derivatives of Eq. (|69ll as 

= ^.Tl^a'^a,)/, (72) 

where the =F sign on the right-hand side is still the inverse of the s-sign. 

We apply these derivatives to the position vector X^ and obtain, as expected, 

— — ±± =u^±. (73 

ds B p ^ ^ 

We consider now the four fractal functions — X'^(s, e^, eg). At this descrip- 
tion level, there is no reason for (— X^)(s, e^, e^) to be everywhere equal to 
— (X'^)(s, e^, £s), owing to a local breaking of the P and T symmetries. We can 
therefore apply, to the two different elementary displacements of these functions, 
the canonical decomposition and the subsequent reasoning applied above to the 
Xf^ function, and obtain the following doublings of the classical velocity: 

^^^^ ~M 

— — ±± =u±±. (74) 
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We use notation analogous to that employed for the case, with the 
addition of a tilde. If we consider the simplest peculiar case when the breaking of 
the symmetry dx^ ^ —dx^ is isotropic as regards the four spacetime coordinates 
(i.e. the signs corresponding to the four ^ indices are chosen equal), we are left 
with eight non-degenerate components - four w± ± and four ± - which can 

be used to define a bi-quaternionic four- velocity. We write 



62 



'++) 



ei 



ez. (75) 



The freedom in the choice of the actual expression for is constrained by 
the following requirements: at the limit when — + dX^ and — > ds, every 
e^-term in Eq. (|75|) goes to zero, and, as = vtj^ in this limit, one recovers 



the complex velocity of Eq. ((T5|) . — 



1+ - - V-+)]/'2-] at the 



classical limit, every term in this equation vanishes, save the real term, and the 
velocity becomes classical, i.e. real: = 

The bi-quaternionic velocity thus defined corresponds to a bi-quaternionic 
derivative operator (f/ds, similarly defined, and yielding, when applied to the 
position vector X^, the corresponding velocity. For instance, the derivative 
operator attached to the velocity in Eq. H75(l is written 




-I 



Substituting Eq. (|72|l and its tilde counterpart into Eq. Il7f)(l . we obtain the 
expression for the bi-quatcrnionic proper-time derivative operator 

— = [1 -f ei + 62 + (1 - i)e:i]— + V^d^ + iVd^'d^, (77) 

the -I- sign in front of the Dalembertian proceeding from the choice of the metric 
signature (-1-,—,—,—). We keep here, for generality, the d/ds term, stressing 
that it is actually of no use, since the various physical functions are not ex- 
plicitly dependent on s. Therefore, the form that we finally obtain for the 
bi-quaternionic derivative operator is unchanged with respect to the previous 
complex operator: this is another manifestation of the covariance of the scale 
relativity description. It is easy to check that this operator, applied to the 
position vector , gives back the bi-quaternionic velocity of Eq. H75() . 

The expression that we have written for in Eq. (|75|l is one among the 
various choices we could have retained to define the bi-quaternionic velocity. 
The main constraint limiting this choice is the recovery of the complex and real 
velocities at the non-relativistic-motion and classical limits. We also choose 
such as to obtain the third term on the right-hand side of Eq. 177|l under the form 
of a purely imaginary Dalembertian, which allows us to recover an integrable 
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equation of motion. To any bi-quatcrnionic velocity satisfying both prescriptions 
there corresponds a bi-quaternionic derivative operator (f/ds, similarly defined, 
and yielding this velocity when applied to the position vector But, whatever 
the definition retained, the derivative operator keeps the same form in terms of 
the bi-quaternionic velocity V, as given by Eq. (|83|l . Therefore, the different 
choices allowed for its definition merely correspond to different mathematical 
representations leading to the same physical result. 



4.2 Bi-quaternionic geodesic equation 

We now apply the generalized equivalence principle, as developed in Sec. 13.21 
The free-motion equation issued from this principle is the geodesic equation 

where is the bi-quaternionic four- velocity, e.g., as defined in Eq. H75|l . 
We introduce the bi-quaternionic action according to 

5S = d^S rfx^ = -mc (5a;^. (79) 

We thus obtain the bi-quaternionic four-momentum, as 

= mc - -9^5. (80) 

We can now define a bi-quaternionic wavefunction, which is again a re- 
expression of the action and which we write 

i^-^d^i: = ^d^s, (81) 

coo 

using, on the left-hand side, the quaternionic product. This gives for the bi- 
quaternionic four- velocity, as derived from Eq. (|80|l . 

V,, = i^i^-^d^^. (82) 
m 

We could choose, for the definition of the wavefunction in Eq. (|81|l . a com- 
mutated expression on the left-hand side, i.e. (9^V')V'~^ instead oiip^^d^ijj. But 
with this reversed choice, owing to the non-commutativity of the quaternionic 
product, we could not obtain the motion equation as a vanishing four-gradient, 
as in Eq. H89|l . Therefore, we retain the above choice as the simplest one, i.e. 
yielding an equation which can be integrated. 



4.3 Free-particle bi-quaternionic Klein-Gordon equation 

We can now replace the covariant derivative operator cf/ds by its expression in 
Eq. H77|l . As we only consider s-stationary functions, i.e. functions which do 
not explicitly depend on the proper time s, this operator reduces to 

i- = V''d,, + iVd''d^. (83) 
as 
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The equation of motion (|78|l is thus written 

{V'd, + iVd^d,) = 0. (84) 
Replacing V^, respectively V, by their expressions in Eq. (|82|l . we obtain 

( i^^-^d^ijd^ + Wd'^d^] U'^dM = 0. (85) 
m \ m J ^ ' 

As in Sec. 12.51 the choice Sq ~ ^rnD allows us to simplify this equation as 

1^-19"^ 9,(1^-19^^) + ^a"9.(V'~'a^V') = 0. (86) 
The definition of the inverse of a quaternion 

V^V^ = = 1, (87) 

implies that '0 and tA^^ commute. But this is not necessarily the case for and 
dfj_ip~^ nor for -0"^ and dfj,ip and their contravariant counterparts. However, 
when we differentiate Eq. (|87|l with respect to the coordinates, we obtain 

^d^i^-^ = -(a^v)V'"' 

^p-'d^^ = -(a^^-i)V, (88) 

and identical formulae for the contravariant analogues. 

Developing Eq. using Eqs. (|^ and the property d^d^dfj, — dfj,d'^d^, we 
obtain, after some calculations, 

a^[(9''a,V)^-i] = 0. (89) 

We integrate this four-gradient and write 

(9"a,V)V'"' +C = 0, (90) 

of which we take the right product by ^ to obtain 

d^d^ij + Cip^Q. (91) 

We therefore recognize the Klein-Gordon equation for a free particle with a 
mass m (so that rn?c^ /h? = C), but now generalized to complex quaternions. 

5 Dirac equation 

We here use a long-known property of the quaternionic formalism, which allows 
to obtain the Dirac equation for a free particle as a mere square root of the 
Klein-Gordon operator (see, e.g., [211 EZ])- We first develop the Klein-Gordon 
equation, as 

1 9^1/) d'^ip d'^tp 9^-0 m'^c^ 

Thanks to the property of the quaternionic and complex imaginary units 
e\ — e\ — e\ — P = —1, we can write Eq. H92|) under the form 



-2^ = ^3^^2 + ^el^^ + e,—e, + ^ -^e,^e,. (93) 
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We see that Eq. H93|) is obtained by applying twice to the bi-quaternionic 
wavefunction ip the operator d/cdt written as 



ld_ 

cdi 



63^62 
ox 



d 
oy 



az 



rac , . 



(94) 



The three Conway matrices )e2, ei( )i and e3( )ei [221, figuring in 
the right-hand side of Eq. (|94|l . can be written in the compact form —a^, with 

cjk 
Ok 

the tTfc being the three Pauh matrices, while the Conway matrix 



63 ( )e3 



is recognized as the Dirac /3 matrix. We can therefore write Eq. (|94|) as the 
non-covariant Dirac equation for a free fermion 














\ 




1 














-1 







V 








-1 





c~dt 



. mc 



(95) 



The covariant form, in the Dirac representation, can be recovered by apply- 
ing 163 ( )e3 to Eq. 

The isomorphism which can be established between the quaternionic and 
spinorial algebrae through the multiplication rules applying to the Pauli spin 
matrices allows us to identify the wavefunction ip with a Dirac spinor. Spinors 
and quaternions are both a representation of the SL(2,C) group ( see Ref. 
for a detailed discussion of the spinorial properties of bi-quaternions) . 



6 Pauli equation 

Finally it is easy to derive the Pauli equation, since it is known that it can 
be obtained as a non-relativistic-motion approximation of the Dirac equation 
[62\ . Two of the components of the Dirac bi-spinor become negligible when 
V << c, so that they become Pauli spinors (i.e. in our representation the bi- 
quaternions are reduced to quaternions) and the Dirac equation is transformed 
into a Schrodinger equation for these spinors with a magnetic dipole additional 
term. Such an equation is just the Pauli equation. Therefore, the Pauli equation 
is understood in the scale-relativistic framework as a manifestation of the frac- 
tality of space (but not time), while the symmetry breaking of space differential 
elements is nevertheless at work. 



7 Conclusion 

Four fundamental motion equations of standard microphysics have been recov- 
ered, in the framework of Galilean scale relativity, as geodesic equations in a 
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fractal space (Schrodinger and Pauli), then in a fractal spacetime (free Klein- 
Gordon and Dirac). It is interesting to note how the change from classical to 
quantum non-relativistic motion, then from quantum non-relativistic to quan- 
tum relativistic motion arises from successive symmetry breakings in the fractal 
geodesic picture. 

First, the complex nature of the wavefunction is the result of the differential 
(proper) time symmetry breaking, which is the simplest effect arising from the 
fractal structure of space (spacetime). It allows a statistical interpretation of 
this wavefunction in the form of the Born postulate (which is demonstrated in 
this framework in the stationary one-dimensional case). At this stage, Galilean 
scale relativity with a complex wavefunction permits the recovery of both the 
Schrodinger (Sec. ^ and Klein-Gordon (Sec. 13 equations. 

To go on with the description of the elementary properties encountered in 
the microphysical world, we consider further breakings of fundamental symme- 
tries, namely the differential coordinate symmetry (dx^ 'r-^ ~dx^) breaking and 
the parity and time reversal symmetry breakings. They provide a four-complex- 
component wavefunction (i.e. a eight component wavefunction), of which the 
most natural mathematical representation is in terms of bi-quaternionic num- 
bers. We therefore obtain the spinorial and the particle-anti-particle nature 
of elementary objects which we can describe as Dirac spinors. The Klein- 
Gordon equation can be recovered from a generalized equivalence principle in 
a bi-quaternionic form which naturally yields the free Dirac equation (Sees. 01 
andEJ. 

The Pauli equation, which is a non-relativistic motion approximation of the 
Dirac equation, then follows. It is understood, in the scale-relativistic frame- 
work, as a manifestation of the fractal! ty of space (not time), while the symmetry 
breaking of the space differential elements is nevertheless at work (Sec. EJ. 

It is worth stressing here that a non-differentiable and fractal spacetime is es- 
sentially non-local, since the "particles" are identified with bundles of geodesies. 
Therefore, we recover the non-locality of the wavefunction of standard quantum 
mechanics. Moreover, having now derived the Dirac equation and the bi-spinor 
nature of the wavefunction exactly in its standard quantum mechanical form 
(there is no missing or additionnal variable in the non-diffcrcntiable spacetime 
representation, but only a change of variables with the same number of de- 
grees of freedom), some profound aspects of quantum mechanics such as the 
EPR paradox and the breaking of Bell inequalities are also recovered in the new 
framework. 

Aknowledgments. The authors thank Dr. T. Lehner and Pr. J. P. Connerade 
for their encouraging support and E. Serie and the referees for useful remarks. 

A Complex numbers and bi-quaternions as a co- 
variant representation 

A.l Introduction 

The first axioms of quantum mechanics state that one defines a state function 
(or probability amplitude) V which is complex, is calculated as "0 = + '02 
for two alternative channels, and is such that the probability density is given 
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by i/'V'^ |34| . Is it possible to derive them from more fundamental principles? 
Such an understanding, impossible in the framework of quantum mechanics 
itself (since these statements are its basic axioms), must be looked for in an 
enlarged paradigm. In scale relativity, we extend the founding stones of physics 
by giving up the hypothesis of spacetime differentiability. One of the main and 
simplest consequences is that the velocity field becomes two-valued, implying a 
two-valuedness of the Lagrange function and therefore of the action. Finally, 
the wavefunction is defined as a re-expression of the action, so that it will also 
be two- valued in the simplest case (see main text). 

But we have, up to now, admitted without justification that this two-valuedness 
is to be described in terms of complex numbers. Indeed, our equations are not 
simply a "pasting" of real and imaginary equations, but involve the complex 
product from the very beginning of the calculations. In particular, the geodesic 
equation for a free particle takes the form of the equation of inertial motion 

^.^^ = 0, (96) 
in terms of the "covariant" complex time derivative operator 

^4 + V,V-,-I,A^ (97) 

Equation amounts to Schrodinger's equation (see main text). Since it 
corresponds to a second derivative, the complex product has mixed the real and 
imaginary quantities in a very specific way, so that it is not at all trivial that 
Schrodinger's equation is recovered at the end. 

The aim of this appendix is to address the question: what is the justification 
of using complex numbers, then bi-quaternions for representing the successive 
doublings of variables? As we shall see, complex numbers achieve a particu- 
lar representation of quantum mechanics in terms of which the fundamental 
equations take their simplest form. Other choices for the representation of the 
two-valuedness and for the new product are possible, but these choices would 
give to the equations a more complicated form, involving additional terms (al- 
though the physical meaning would be unchanged). Bi-quaternions follow as a 
further splitting, at another level, of the complex numbers. 



A. 2 Origin of complex numbers and quaternions in quan- 
tum mechanics 

Let us return to the step of our demonstration where complex numbers are 
introduced. All we know is that each component of the velocity now takes two 
values instead of one. This means that each component of the velocity becomes a 
vector in a two-dimensional space, or, in other words, that the velocity becomes 
a two-index tensor. So let us introduce generalized velocities 

^(V'',U'') ; fc = l,2,3 ; a = -,+. (98) 

This can be generalized to other physical quantities affected by the two- 
valuedness: scalars A of the position space become vectors Aa- of the new 2D- 
space, etc. While the generalization of the sum of these quantities is straigh- 
forward, = A^ + B^, the generalization of the product is an open question 
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at this stage. The problem amounts to finding a generahzation of the standard 
product that keeps its most fundamental physical properties (e.g., internal com- 
position law), when some of them cannot escape to be lost (e.g., commutativity 
when jumping to quaternions). 

From the mathematical point of view, we are here confronted with the well- 
known problem of the doubling of algebra (see, e.g., 23 )■ The effect of the 
symmetry breaking dt ^ —dt (or ds —ds) is to replace the algebra A, 
in which the classical physical quantities are defined, by a direct sum of two 
exemplars of A, i.e. the space of the pairs (a, 6) where a and b belong to 
A. The new vectorial space A^ must be supplied with a product in order to 
become itself an algebra (of doubled dimension). The same problem shows 
again when one takes also into account the symmetry breakings dx'^ <-*■ —dx^ 
and x^ ^ ~x^: this leads to new algebra doublings. The mathematical solution 
to this problem is well known: the standard algebra doubling amounts to supply 
A^ with the complex product. Then the doubling of IR is the algebra (D of 
complex numbers, the doubling (D^ of (D is the algebra H of quaternions and 
the doubling of quaternions is the algebra of Graves-Cayley octonions. The 
problem is that the iterative doubling leads to a progressive deterioration of the 
algebraic properties. Namely the quaternion algebra is non-commutative, while 
the octonion algebra is also non-associative. But an important positive result 
for physical applications is that the doubling of a metric algebra is a metric 
algebra 

These mathematical theorems lead us to use complex numbers, then quater- 
nions, in order to describe the successive doublings due to discrete symmetry 
breakings at the infinitesimal level, which are themselves more and more pro- 
found consequences of spacetime non-differentiability. However, we give in what 
follows complementary arguments of a physical nature, which show that the use 
of the complex product in the first algebra doubling has a simplifying and co- 
variant effect. 

In order to simplify the argument, let us consider the generalization of scalar 
quantities, for which the product law is the standard product in JR. 

The first constraint is that the new product must remain an internal com- 
position law. We also make the simplifying assumption that it remains linear in 
terms of each of the components of the two quantities to be multiplied. There- 
fore, we consider a general form for a bilinear internal product 

C = A" (99) 

where the matrix 17^^^ is a tensor (similar to the structure constants of a Lie 
group) that defines completely the new product. 

The second physical constraint is that we recover the classical variables and 
product at the classical limit. The mathematical equivalent is the requirement 
that A still be a sub-algebra of A^ . Therefore, we identify E A with (aq, 0) 
and we set (0, 1) = a. This allows us to write the new two-dimensional vectors 
in the simplified form a = a^} + aia, so that the product is now written 

c ~ (oo + aio) {bo + bia) = ao&o + aibio^ + [a^bi + aibo)oi. (100) 

The problem is now reduced to finding a^, i.e. only two Q. coefficients instead 
of eight. 

= + ijJia. (101) 
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Let us now return to the beginning of our construction. We have introduced 
two elementary displacements, each of them consisting of two terms, a classical 
part and a fluctuation (see Eq. (|ll|l ) 

dX+{t) = v+ dt + d^+{t), 

dX_{t)=v^ dt + dS,-{t). (102) 

Therefore, one can define velocity fluctuations w+ = d£,+/dt and w_ = 
d£,-/dt, then a complete velocity in the doubled algebra [201 

V + W = a + a . (103) 



In terms of standard methods, this writing would be forbidden since the 
velocity W is infinite from the viewpoint of usual differential calculus (it is 
cx dt~^/'^). But we give meaning to this concept by considering it as an explicit 
function of the differential element dt, which becomes itself a variable. 

Now, from the covariance principle, the Lagrange function in the Newtonian 
case should strictly be written: 

£ = im Cl{{V + Wf) = im {Ci{V^) + Ci{W^)) (104) 

We have C£{W) — 0, by definition, and Cl{VW) — 0, because they are mu- 
tually independent. But what about 'Cl{W'^) ? The presence of this term would 
greatly complicate all the subsequent developments toward the Schrodinger 
equation, since it would imply a fundamental divergence of non-relativistic quan- 
tum mechanics. Let us expand it 

4C?(>V2) = 'Cl{[{w++w^)-a{w+-w^)Y) 

= Cl{{w\ + w'i){l + o?) - 2a{wl - w^_) 

+ 2w+w_(l - a^)). (105) 

Since Ci{w\) = CKvPl) and C£{w+w^) — (they are mutually indepen- 
dent), we finally find that C£{W'^) can only vanish provided 

= -1, (106) 

namely a — ±i, the imaginary. Therefore, we see that the choice of the complex 
product in the algebra doubling plays an essential physical role, since it allows 
the suppression of additional infinite terms in the final equations of motion. 



A. 3 Origin of bi-quaternions 

A last point to be justified is the use of complex quaternions (bi-quaternions) 
for describing the new algebra doublings that lead us to bi-spinors and the Dirac 
equation. One could think that the argument given in Sec. lA. 21 implies the use 
of Graves-Cayley octonions (and therefore the giving up of associativity) in the 
case of three successive doublings as considered in this paper. However, these 
three doublings are not on the same footing from a physical pointof view: 

(i) The first two-valuedness comes from a discrete symmetry breaking at the 
level of the differential invariant, namely dt in the case of a fractal space (yielding 
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the Schrodinger equation) and ds in the case of a fractal spacetime (yielding the 
Klein-Gordon equation). This means that it has an effect on the total derivatives 
d/dt and d/ds. This two-valuedness is achieved by the introduction of complex 
variables. 

(ii) The second two-valuedness (differential parity and time reversal vio- 
lation, "dX"), which is specifically introduced and studied in the present pa- 
per, comes from a new discrete symmetry breaking (expected from the giv- 
ing up of the differentiability hypothesis) on the spacetime differential element 
dx^^ ^ —dx'^. It is subsequent to the first two-valuedness, since it has an effect 
on the partial derivative = d/dx^ that intervenes in the complex covariant 
derivative operator, namely, 



(iii) The third two-valuedness is a standard effect of parity (P) and time reversal 
(T) in the relativistic-motion situation, which is not specific to our approach 
and is already used in the standard construction of Dirac spinors [S^]. It does 
not lead to a real information doubling, since Dirac spinors still have only four 
degrees of freedom as Pauli spinors do. Therefore, from the second and third 
doublings, complex numbers, then quaternions, can be introduced, which will 
affect variables which are already complex due to the first, more fundamental, 
doubling. This leads to the bi-quaternionic tool we use here. 

It is worth noting that these symmetry breakings are effective only at the 
level of the underlying description of elementary displacements (namely in the 
non-diffcrcntiable fractal spacetime) . The effect of introducing a two-valuedness 
of variables in terms of double symmetrical processes amounts to recover sym- 
metry in terms of the bi-process, and therefore in terms of the quantum tools 
which are built from it. In reverse, it opens a possible way of future investigation 
into the origin of other features characteristic of the microphysical world. 

A. 4 Complex representation and covariant Euler-Lagrange 
equations. 

Let us confirm by another argument that our choice of a representation in 
terms of complex numbers for the two-valuedness of the velocity (namely in 
the minimal situation that leads to the Schrodinger equation) is a simplifying 
and covariant choice. We demonstrate hereafter that the standard form of the 
Euler-Lagrange equations is conserved when we combine the two velocities in 
terms of a unique complex velocity. 

In a general way, the Lagrange function is expected to be a function of the 
variables x and their time derivatives x. We have found that the number of 
velocity components x is doubled, so that we could have written 



(107) 



L = L{x, X-, t). 



(108) 



Instead, we have made the choice to write the Lagrange function as 



L = L{x,V,t). 



(109) 



25 



We now justify this choice by the covariance principle. Re-expressed in terms 
of x+ and i_, the Lagrange function is written 



L = L \ X, — - — x+ H — a;_ , t 



Therefore we obtain 
dL 



1-i dL 
dV 



dL 



1 + i dL 
9V' 



di+ 2 dV ' ai_ 2 
while the new covariant time derivative operator reads 

1 



dt 



i d 
~ dtl 



1 + i d 
2 dtl 



(110) 



(111) 



(112) 



Let us write the stationary action principle in terms of the Lagrange function 
ofEq. fTUSI) 



SS^S 



L{x,x+,X-,t) dt = 0. 



It becomes 



t2 



dL 



5x 



dL 



5x4 



dL 



5x- 



dt = 0. 



dx dx+ dx- 
Since 5i+ = d{5x)/dt^ and 6x- = d{Sx) / dt^ , it takes the form 



— 6x — 
dx dV 



1-i d 
2 dtl 



1 + i d 
2 dH 



Sx] dt^ 0, 



I.e. 



g^Sx + --Sx]dt = 0. 



(113) 



(114) 



(115) 



(116) 



The subsequent demonstration of the Lagrange equations from the stationary 
action principle relies on an integration by parts. This integration by parts 
cannot be performed in the usual way without a specific analysis, because it 
involves the new covariant derivative. 

The first point to consider is the Leibniz rule for the covariant derivative 
operator d /dt. Since d /dt = d/dt + V . V — iDA is a linear combination of first- 
and second-order derivatives, the same is true of its Leibniz rule. This implies 
an additional term in the expression for the derivative of a product (20j : 



d fdL 



dt \dV 



6x 



d dL ^ dL d ^ ^.^„dL 

— — .6x + — . —Sx - 2tW— .VSx. 

dt dV dV dt dV 



(117) 



Since Sx{t) is not a function of x, the additional term vanishes. Therefore, 
the above integral becomes 



t2 



dx ~ dt dv + dt [dv 



dt = 0. 



(118) 
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The second point is the integration of the covariant derivative. We define a 
new integral as being the inverse operation of the covariant derivation, i.e. 



j <^f^f (119) 



in terms of which one obtains 



t2 



- 0, (120) 



since Sx{ti) — Sx(t2) = by definition of the variation principle. Therefore, the 
action integral becomes 

\dx dt dV ) ^ ' 

And finally we obtain generalized Euler-Lagrange equations that read 

= (122) 
dt dV dx ^ ' 

They take the form obtained by writing a stationary action principle based 
on Eq. 1109|) . Moreover, once the transformation d/dt d /dt is done, this 
form is nothing but the standard classical form. This result reinforces the identi- 
fication of our tool with a "quantum-covariant" representation, since this Euler- 
Lagrange equation can be integrated in the form of a Schrodinger equation. 
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